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Abstract. Utrecht has a long tradition in both spectroscopy and mass-loss studies. 
Here we present a novel methodology to calibrate mass-loss rates on purely spectro- 
scopic grounds. We utilize this to predict the final fates of massive stars, involving 
pair-instability and long gamma-ray bursts (GRBs) at low metallicity Z. 



1. Introduction 

Mass loss is an important ingredient in massive star evolution modelling. Nowhere is it 
more dominant than for the most massive stars, which are thought to evolve chemically 
homogeneously (Grafener et al. 2011). Mass loss determines the final fate. This may 
involve normal supernovae (SNe) or pair-instability SNe, leaving either normal black 
holes or intermediate mass-black holes (IMBHs) behind (Belkus et al. 2007; Yungelson 
et al. 2008). The point is that very massive stars (VMS) up to 300M o - and perhaps 
even as high as ~1000M o through collisions in clusters - are now thought to exist in 
nature (Crowther et al. 2010; Bestenlehner et al. 2011). 

In recent years a debate has arisen regarding the roles of stellar wind versus erup- 
tive mass loss, as winds have been found to be clumped, which resulted in the reduction 
of empirical mass-loss rates. Stellar evolution modellers however generally employ 
theoretical mass-loss rates. These are already reduced by a moderate factor (of ~2-3) 
compared to non-corrected empirical rates. A key question is whether these reduced 
rates are correct or if they need to be reduced even more. 

Whilst stationary winds in O and Wolf-Rayet (WR) stars are ubiquitous, it is not 
clear if objects like tj Car have experienced a special evolution, e.g. involving a merger, 
or if a majority of massive stars would encounter eruptive mass-loss phases. Alterna- 
tively, for the most massive main-sequence WNh stars there is strong evidence for a 
Eddington parameter T-dependent mass loss (Grafener et al. 2011). In other words, for 
VMS the role of stationary mass loss has increased in recent years. 

We recently introduced the transition mass-loss rate M trans between O and WR 
stars (Vink & Grafener 2012). Its novelty is that it is model independent. All that is 
needed is the spectroscopic transition point in a data-set, and to determine the stellar 
luminosity. This is far less model-dependent than conventional mass-loss diagnostics. 
Our results suggest that the rates provided by Vink et al. (2000) are of the right order 
of magnitude in the 50M G range, but alternative mechanisms might be needed at lower 
Z, particularly relevant for the occurrence of GRBs. 
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Figure 1 . Mass-loss predictions versus the Eddington parameter F - divided by 
iVr 1 . Symbols correspond to models of different mass ranges (Vink et al. 201 la). 



2. The mass loss versus Eddington-Gamma dependence 

Vink et al. (2011a) discovered a kink in the slope of the mass-loss vs. Y relation at 
the transition from optically thin O-type to optically thick WR-type winds. FigureQ] 
depicts mass-loss predictions for VMS as a function of the Eddington parameter T. For 
ordinary O stars with "low" Y the M x Y x relationship is shallow, with x —2. There 
is a steepening at higher Y, where x becomes -5. Here the optical depths and wind 
efficiencies exceed unity. This result from Monte Carlo modelling, i.e. that the O to 
WR mass loss transition point coincides with n - 1, can also be found analytically. 



2.1. The transition mass-loss rate 

Lamers & Cassinelli (1999) provided momentum considerations for dust-driven winds 
that we apply to line-driven winds. The integral form of the momentum equation has 
four terms, but as hydrostatic equilibrium is a good approximation for the inner wind, 
and the gas pressure gradient is small beyond the sonic point, the 2nd and 3rd terms are 
neglected, leading to: 

f°° , dv r°° GM , 

Anr 2 pv—dr+ —-{\-Y)pAnr 2 dr = 0. (1) 

JR* dr Jr s r 2 

When we utilize the mass-continuity equation M = 4nr 2 pv, we retrieve 

M—dr - M Voo = 4nGM (Y(r) - \)pdr (2) 

where r s represents the sonic radius and Y(r) = is the Eddington factor with 

XCO 
/epp dr, 

we obtain 
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If we assume that Y is larger than unity, the factor (T - l)/r is close to one, which leads 
to Mvoo = L/ct. We can now derive the unique wind efficiency condition n - ^jjr - 
t - 1 for the transition point from optically thin O-star winds to optically-thick WR 
winds, and we retrieve a model-independent M. In case we have an empirical data-set 
available that includes luminosity determinations we can provide the transition mass- 
loss rate M tra ns> simply from the transition luminosity L trans (and the terminal velocity 
Voo): Mtrans = . This transition point is obtainable just by spectroscopic means, and 
independent of any wind-clumping assumptions. 

2.2. Do very massive stars make IMBHs or pair instability SNe? 

To address these questions, we require accurate mass-loss rates for VMS. As Vink & 
Grafener (2012) were able to calibrate the Vink et al. (2000) mass-loss prescriptions 
at the high-mass end, and given that these rates agree well with the Crowther et al. 
(2010) rates for the 30 Dor WNh core stars for moderate clumping factors D ~ 10, we 
employ the Vink et al. (2000) rates for now. When starting with 300M o , we find M 
= 10~ 4 - 2 M G yr _1 . For a 2.5Myrs lifetime this leads to a total main-sequence mass lost 
of -150M G . Extra wind mass loss during the core helium phase will further evaporate 
these stars. There seems little space left for eruptive mass loss. Our results would imply 
that IMBHs and pair-instability SNe are unlikely, unless we move to lower Z. 

3. GRBs and the metallicity dependence of Wolf-Rayet stars 

Mass loss at low Z has gained attention due to the issue of cosmic reionization by pop- 
ulation m stars. Also, massive stars are thought to be the progenitors of GRBs. Within 
the collapsar model, GRB progenitors require 2 properties: (i) rapidly rotating cores, 
and (ii) the absence of hydrogen envelopes. Therefore, GRB progenitors are thought 
to be rotating WR stars. A potential pitfall is that WR star have high mass loss which 
could remove all the angular momentum before core collapse. In the rapidly rotating 
stellar models of Yoon & Langer (2005), the objects evolve "quasi-homogeneously". 
The objects are subjected to a strong magnetic coupling between the stellar core and 
envelope. If rapid rotation can be maintained as a result of lower main-sequence mass 
loss in lower Z galaxies, the objects may avoid spin-down during a RSG or LBV phase, 
and directly become rotating WR stars. If WR winds also depend on Fe driving (Vink 
& de Koter 2005), the objects could remain rapid rotators towards the end of their lives, 
enabling GRB formation, but exclusively at low Z. 

Recent GRB observations however suggest that GRBs are not restricted to low Z. 
Indeed, there seems to be a need for a high Z GRB channel. We have recently identified 
a sub-group of rotating Galactic WR stars. This would allow for a potential solution 
to the problem (Vink et al. 2011b; Grafener et al. 2012). Spectropolarimetry surveys 
show that the majority of WR stars have spherically symmetric winds indicative of slow 
rotation, but a small minority display signatures of a spinning stellar surface (see Fig. 2 
for the example of WR 134). We recently found the spinning subgroup surrounded by 
ejecta nebulae - thought to be ejected during a recent RSG/LBV phase - suggesting 
that these objects are still young, and rotating. 
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Figure 2. Galactic Wolf-Ray et star WR 134. (a) Position Angle (PA) of the po- 
larisation, (b) Degree of linear polarisation, (c) Stokes I "intensity". 

If core-surface coupling were strong enough the cores would not rotate rapidly 
enough to make GRBs. However, if core-envelope coupling is less efficient stars may 
have sufficient angular momentum in the core to make a GRB. At high Z, the objects 
would in most cases still spin down due to mass loss, but with our post-RSG/LBV sce- 
nario one would not exclude the option of high Z GRBs. Still, low Z remains preferred 
due to weaker WR winds (Vink 2007). 
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